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1 Introduction 



The setting of group field theory (GFT) produces tensor models which generalize matrix 
models to structures of dimensions higher than two. In particular, the Boulatov model pQ 
generates a sum over topologies of 3-manifolds (as well as other objects) weighted by the 
corresponding Ponzano-Regge amplitude, which is a model for the partition function of 3d 
Riemannian quantum gravity. As a field theory, it provides a good framework to address the 
renormalization of theories where one has to sum over topologies of the spacetime manifold 
(or cell complex). In particular, it stands as a promising approach to quantum gravity, via 
spin foam models, where one aims at a background free quantization of general relativity, 

Several interesting results and advances in GFTs have been obtained during the last 
years. From contacts with simplicial geometry [6] to mechanisms for emergent matter fields 
(incorporating a new/extra gauge symmetry) [7], from contacts to noncommutative geometry 
[8] and quantum groups [9], [10] to classical geometrodynamics extraction and Bose condensate 
techniques pj] , indeed GFT evolves fast and in several directions. 

From the quantum field theory perspective, efforts have been concentrated so far mostly 
on power-counting results, either in the Boulatov model and its higher dimensional extensions 
(coined flat spin foam model) p2|-[T7], or in theories hopefully related to quantum gravity, 
[IH[T9]. An exciting development of these power-counting results is the recently defined large 
N limit, [20], for colored GFTs [21]- [23] where the analog of the planar graphs of matrix 
models are found to be 2-complexes corresponding to the 3-sphere. 

Still, one may be interested in making sense of the theory without the 1/N expansion. A 
key feature which remains to be unraveled is a locality principle - typically in ordinary scalar 
field theory, it says that renormalization appears when probing high internal scales, i.e. short 
distances, with low external momenta. For GFTs, a generalized locality principle was seen 
to emerge in [18]. Also in [16], the flat spin foam model was considered as a generalization 
of the zero coupling limit of 2d Yang-Mills theory, [25], to cell 2-complexes. Thus, the set 
of interest which was used to localize the integrals was the set of flat connections on the 
complex, and the divergence degree could be extracted in terms of a discrete analog of the 
twisted De Rham co homology for the covariant differential (1a, A being a flat connection. 

In this paper, we will use the setting which proved fruitful in [IB] to compute radiative 
corrections to the propagator. It is worth emphasizing the main observation which makes 
this method efficient. The manifold in the model is (several copies of) a compact Lie group, 
say SU(2). Thus, the momentum space is described in terms of half-integer spins (irreducible 
representations of SU(2)). To understand the localization in the Boulatov model, it is much 
simpler to explicitly perform the sums over the spins, which produce delta distributions 
on the manifold. Getting a similar picture of the localization without the explicit result of 
summing over the spins is certainly a challenging task, and it may help to extend the analysis 
to other models (where the integrals may be concentrated around something else than flat 
connections) . 

To make our program on the 2-point function clear, let us remember the kind of expansion 
usually performed in the ordinary <p\ theory, and which leads to the key ideas of mass 
subtraction and wave function renormalization. The standard textbook procedure is to re- 
express the bare quantities in terms of the physical ones, say the field: 0& = Z l / 2 <p r , and the 
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mass m 2 = m 2 Z — 5 m , so that the Lagrangian, written with renormalized quantities, picks 
up counter-terms of the form: 

P 2 5z-Sm, (1) 

in momentum space, with 5z — Z — 1. Obviously, the term proportional to p 2 comes from 
the part of the action with derivatives of the field, 0<9 2 in direct space. Then, 5 m ,5z 
are extracted from radiative corrections to the 2-point function, using some renormalization 
prescription. This is done by evaluating the graph amplitude at a given external scale, like 
Pcxt = m ri f° r 3m, and evaluating its derivatives with respect to p 2 at this scale for 5z- 

Note that in [19] the mass subtraction has been roughly worked out for the EPRL/FK 
model, a candidate for quantum gravity via spin foams. Nevertheless, we would like to under- 
stand first how the prime concepts of leading divergences and counterterms find themselves 
a natural extension in the simpler Boulatov-Ooguri tensor model. 

This model differs from the usual 4>\ theory by the fact that the quadratic part of the 
Lagrangian does not contain derivatives of the field, but only a mass term, 

J m 2 2 + A 4 . (2) 

oo, roug the propagator behaves like: C = J dp^2~. Consequently, we will look for 
organizing the expansion of each 2-point graph, not as an asymptotic expansion with the 
cut-off, but as a Taylor expansion where the n-th order is proportional to a n-th derivative 
of the propagator. Symbolically: 

A fe Q + A *i C (l) + A k 2 C {2) + _ _ _ ^ (3) 

where A is the ultraviolet cut-off. We are interested only in the leading orders of this 
expansion, i.e. the terms which diverge when the cut-off A goes to infinity. A divergence 
with the cut-off at the zeroth order, i.e. ko > 0, is a constant term in momentum space 
(given the specific form of the propagator), and gives a mass renormalization 5 m . We expect 
that there is no term with by symmetry. The term proportional to behaves like p 2 
in momentum space. Thus if the corresponding leading term diverges with A (i.e. k 2 > 0), 
we should get a wave-function renormalization 5z- 

This is indeed what we will find, from a graph at the second order in the coupling constant. 
However, due to the form of the initial action (J2]), there is no possible counter-terms for the 
divergences with C {2) . This is the sign that the action needs an addition of a quadratic term 
with second derivatives of the field, so that it takes the form: 

J 0A0 + m 2 2 + A0 4 , (4) 

where A is the Laplace operator on the group manifold. It turns out that the resulting group 
field theory has already been considered in the literature by Oriti ([7j and more references 
therein) with different motivations. 

The organization of the paper is as follows. In Section |2l we define the Boulatov -Ooguri 
tensor model. 

1 The peculiarities of the GFT (being non-local, the covariancc projecting the field onto a rotationally 
invariant sector) play no role in this discussion. 
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• In Section [31 we study a correction to the 2-point function, at the second order in 
the coupling constant, and observe two leading divergences. We show that there is 
indeed a mass renormalization. We also find a divergence proportional to the second 
derivatives of the propagator, which cannot be balanced with a counter-term from the 
initial action. The structure of the corresponding term is studied with details in the 
3d case. 

• In Section HI we study generic graphs of the 2-point function. We obtain an expression, 
in Lemma [3], which shows that the integrand is localized on the set of flat connections 
on the 2-complex, so that the generalized Laplace approximation of [IB] applies. The 
first order in particular is an averaging of the propagator with non-trivial insertions on 
flat connections. 

• In Section El we show that an expansion of the form is directly available for a class 
of 2-point graphs which are simply connected 2-complexes. We then find that there 
exist divergences to the order n in the derivatives of the propagator with n arbitrarily 
high, which are computed from Gaussian moments of the zero order. 

• Section [6J discusses the addition of specific quadratic terms with second derivatives of 
the field in the action. We extract the propagator, and argue that we have a natural 
setting for a scale analysis. 

The technical tools we use do not rely on the specific form of the graphs of the GFT, 
and may be applied to some cellular complexes as well. Nevertheless, there is no hope to 
get an expansion like ([3]) for the generic GFT graphs. The kind of graphs we will consider 
throughout the paper contains all graphs of the colored model [21], a subset of the generic 
graphs, and certainly graphs which are not of the GFT type. 



2 The action and the graph amplitudes 

We start by giving the action and some properties of Feynman graphs for the Z)-dimensional 
GFT of Boulatov-Ooguri type, [T], over a compact Lie group G. Typically, G = SU{2) is 
considered throughout the text. Then in the subsequent sections, we proceed to a specific 
expansion of the 2-point function. 

Fields belong to the Hilbert space of square integrable functions on G D , namely H = 
L 2 (G D ). These fields are actually restricted the so-called gauge invariant fields, in the sense 
that they are invariant under the diagonal right action of the group: 

(j)(gih,g 2 h, ...,g D h) = <j>(g u g 2 ,. ..,5-d) , V7j G G . (5) 

The quadratic part of the action has the form: 

D 

/ IJ dg s <f>(gi, go)(f>(gi, ■ ■ ■ , go) , (6) 

J s=l 
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where dg denotes the Haar measure on G. Due to the restriction (JSJ) to a subset of fields, 
the propagator is not just the inverse of this quadratic part. Instead, one has to use a 
(normalized) degenerate Gaussian measure d/j,c[4>], of covariance C defined by: 

C({gi, . . . , g D }\ {h, • • • , 9d}) ■= I dh \\5{g s hg s x ) . (7) 

J 8=1 

The integral on h is a group averaging which enforces the condition (JSJ). 
The interaction is non-local, of degree (D + 1), namely a theory: 

S in t[(j>] := A J Y\dg a b 0i2 ) i3 v .. )1 d+i ^(£1+1)1,(13+1)2 ,...,(D+1) D 4>D D + 1 ,D 1 ,...,D D ~ 1 • • • 



D+1 

l) 3 4 ,3 5 ,...,3 D + 1 ;S 2 023,2 4 ,...,2 D + 1 ,2 1 ] J Kd a b 9\p) 



where the shorthand notation (f)(gia,gib, . . . ) = <f) la ^ is understood. Each integration vari- 
able appears in two copies of the field. 
The partition function is: 



Z(\) ■= J dnc[<j>\ e- XSinM] . (9) 




Figure 1: D dimensional model propagator and vertex. 

The propagator can be drawn with D strands (see Figured]), each of them corresponding 
to a delta function in (J7|). As for the vertex, each of the (D + 1) legs has D strands. One 
strand goes between two legs, and there is a single strand between any pair of legs. The 
Feynman rules produce Feynman graphs which are actually 2-complexes. Indeed, in addition 
of the set of vertices and edges (or lines, due to propagators), one identifies faces as regions 
being bounded by closed strands. An open strand gives an open face. The corresponding 
amplitudes take the form: 

M{g? t h{~g7 t }) = 



n 



n 5 U ext 



/eF cxt (S) 



n k 

leaf 



(9. 



~ext \ — 1 




(10) 
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where the group elements {g , f xt }^ = i and {g° xt }fLi are the arguments on the external legs, 
L(Q) denotes the set of (internal) edges or lines, F ext ((?) and F mt {Q) the sets of external 
(open) and internal (closed) faces, respectively, of the graph Q. The product over u l G df" 
means the product over the set of lines I belonging to the boundary of the face /. The sign 
eif = ±1 is equal to +1 if the orientations of the face and of the line coincide, to —1 if not 
and to if the line does not belong to the boundary of /. Note that we take, as a convention 
and without loss of generality, a particular orientation of the external leg arguments. 

There is a natural geometric interpretation which generalizes the two-dimensional in- 
terpretation of graphs of matrix models to higher-dimensional cell complexes. Indeed, the 
specific pattern of the vertex mimics the structure of a D-simplex. Each leg is considered to 
be one of the (D + 1) boundary (D — l)-simplices. A strand represents a (D — 2)-simplex on 
the boundary, shared by exactly two (D — l)-simplices. In this view, the propagator is seen 
as a (D — l)-simplex. Furthermore, the Feynman rules produce all possible ways of gluing a 
set of .D-simplices along their boundary. 

Since products of distributions are not always well-defined, formulae such as ffTU]) suffer 
from what one should consider from the quantum field theoretic point of view as divergences. 
In the present situation, they will be referred to as "ultraspin" divergences. 

3 First radiative corrections to the 2-point function 

In three dimensions, a graph is built with three strands per propagator and the vertex is of 
the type 4 , see Figure EJ 



«3 



6 3 

si 
si 



Figure 2: 3d model propagator and vertex. 



3.1 The first 2-point graph and its regular izat ion 

We begin our study with a 2- vertex graph of the 2-point functional (see Figure |3]). The 
combinatorial ingredients of the graph are: three internal edges, L(Q) = {h,l2,h}, three 
external faces, F cxt (£) = {/°, /°, / 3 }, and three internal faces, F int (£) = {/ 12 , Ji 3 , / 23 }. 
Explicitly, the amplitude built from the graph formally reads: 



,,3 3 

Ag({g s }; {~g s }) = Udht ]J5 (gshsCgs)- 1 ) J] 6 (h.hj 1 ) 

J 1=1 s=l Ki<j<3 



2 In the ordinary GFT, there are lots of other 2- vertex graphs. But it is worth observing that the graph 
we consider is the simplest one in the colored GFT of [2"T] . 
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Figure 3: 2-point function. 



To make sense of the product of Dirac delta in pip , we change each of them with a heat 
kernels at time t on G. Thanks to the Peter- Weyl theorem, it can be expanded onto SU(2) 
irreducible representations, labelled by half- integers (spins), 

K t {g)= ^(2j + l)e 4C ^xj{9)- (12) 
jen/2 

Here C(j) is the Casimir in the representation of spin j, and Xj the character. When t goes 
to zero, this goes to + l)xj(g) = $(9), which is indeed the usual expansion of 5(g). 

For small times t, the heat kernel is localized around zero, and its behaviour is close to 
that of the Euclidean kernel. In a neighborhood of the identity, 

K t (g) ~ o Af mG e"^ , (13) 

where \g\ is the Riemannian distance from the identity to g, and 

A t = (4 7 rt)- 1 / 2 . (14) 

Results concerning the divergence degree of graphs in this theory can be found in the 
literature either using this heat kernel regularization, or using a sharp cut-off A on the spin 
expansion of the delta function, 5\(g) = J2f=o(2j + l)Xj(fl')- The cut-off A t is chosen so 
that the divergence degree of a given graph is the same if computed using 5\. Hence, in a 
nutshell, high spins correspond to small t. 

As we explain in the paragraph below, we will apply a saddle point approximation to 
evaluate the graphs. However, we will be interested only in some specific terms which 
diverge with the cut-off A f . For our purpose the asymptotic behaviour (1131) will be sufficient. 
However to perform the full asymptotical expansion^ in exponents of A t , one would have to 
look at higher orders in the expansion of the heat kernel. 

3 Thc heat kernel is the solution of: (d t — A) K t (g) — 0, in which A is the Casimir-Laplace operator on 
G, and with the initial condition lim f _,.o K t (g) = 5(g). 

4 From a field theory perspective, such a full asymptotical series is not really relevant, since it will 
strongly depend on the chosen regularization. The most natural thing to do instead is to go within a 
multiscale analysis, in which one gets interested into subgraphs carrying higher scales than the rest of the 
graph. 
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3.2 Extracting the divergences: Mass and wave-function renor- 
malizations 



The regularized amplitude can be expressed as 



A G;t({9s}; {g s }) = / dh x dh 2 dh 



IPOA&r 1 ) Kt^h- 1 ) Kt^h- 1 ) K t {h 2 h- 1 ) . 

s=l 

(15) 

Changing of variables fcj = h 1 1 hi, for i = 2,3, and using the translation invariance of the 
Haar measure allow us to rewrite the same amplitude as 

Ag; t ({g s }; {g s }) = 



dk 2 dk 3 



dhi5 (gihi(gi) x ) JJ<5 {g s hik s {g s ) x ) 

s=2 

K t (k 2 )K t (k 3 )K t (k 2 k^) . (16) 



It is worth emphasizing that this change of variables is a gauge fixing procedure, whose goal 
is to put to the identity one of the h s (here h s = hi) in the internal faces. Indeed, the 
diagonal invariance of the field ([5]) induces an invariance under an action of G at each node 
of each graph. Like in [27], this gauge invariance can be (partially) fixed by contracting a 
maximal tree in the graph. The new point here is that this procedure has to be implemented 
on graphs with external legs. As it can be observed in the above example, the net result 
is that hi indeed disappears from the internal faces, but then appears in the three external 
faces. We will generalize this procedure to any graph of the 2-point function in the next 
paragraph. 

The reader may worry that the symmetry between the three strands is not explicit any- 
more after this change of variables. Still the symmetry will be completely explicit in our 
final result. 

Using the Gaussian approximation ( fl3l) . one gets, 



A?;t({#J; {g s }) - A t 9 J dk 2 dk 3 C({gi, g 2 , g 3 }; {g x , g 2 k 2 , g 3 k 3 }) e « 



(w 



1*3 I 



|fc 2 fc 3 - 1 | 2 ) 



(17) 

Here the integral over hi has been reabsorbed to form the propagator ([7j), with insertion of 
k 2 , k 3 in the strands 2 and 3. 

The saddle point - We are now in position to perform a saddle point approximation 
around: 

k 2 = k 3 = 1 . (18) 

The fact that there is a single saddle point can be understood from a larger perspective, 
as we explain in Sections HI [5j Mainly, the amplitude of a graph in the Boulatov-Ooguri 
model is a special case of the flat spin foam model considered in [16], where it is shown that 
the integral is generically localized around the set of homomorphisms from the fundamental 
group of the 2-complex to the group G. Here, Q is simply connected, so there is a single 
homomorphism to SU(2). 

The saddle point approximation proceeds by expanding k 2 = e X2 into powers of the Lie 



algebra element X 2 , and the same for k 3 



The integrals over k 2l k 3 are changed to 
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integrals on the tangent space, the Lie algebra su(2), equipped with the Lebesgue measure 
d 3 X 2 d 3 X 3 . The Haar measure is written as: 

dk 2 dk 3 = MX 2 ,X 3 ) d 3 X 2 d 3 X 3 . (19) 

The Hessian - Similarly, the function (|fc2| 2 + \k 3 \ 2 + \k 2 k 3 1 \ 2 ) in the exponential is 
expanded to extract the Hessian, 

N 2 + N 2 + IM^ T = x 2 + x\ + (x 2 - x 3 f + S> 3 (X 2 , X 3 ) . (20) 

The squared quantities in this formula correspond to the norm of the Killing form (in our 
simple SU(2) case, this is the Euclidean norm on R 3 ). We have packed all terms of order 3 and 
higher into the remainder S> 3 (X 2 , X 3 ). Quite clearly, the quadratic part is non-degenerate, 
its kernel being {X 2 = X 3 = 0} . 

In an orthonormal basis for (X 2 ,X^) the Hessian matrix is: 

Hess=f 2 ( 3 ~{ 3 Y (21) 



-h 2J 3/ 

where I 3 is the 3x3 identity matrix. Up to a constant, its inverse is: 

Hess -1 = const x ( 2 / 3 / 3 J . (22) 
\ h 2 1 3 J 

Expanding the propagator - Since it is a product of Dirac delta, we first consider 
derivatives of delta. To avoid test functions, we will regard them through the Fourier expan- 
sion. Remember that the character Xj hi the representation of spin j is the trace Tr^ on the 
carrier space of dimension (2j + 1). This provides us with a matricial picture, in which we 
can write, at least formally: 

00 1 

^ e *) = E - { E w + x ) Tr ^ (9 w n ) • ( 23 ) 

n=0 ' jGN/2 

Next, pick up any basis (Jj)i=i,2,3 of the Lie algebra (anti-hermitian generators), and setEl 
X = X 1 J u so that ( 1231) becomes 

^22, Y ix X i2 ■ ■ ■ X in 

5{ge x ) = ~< E( 2 ^' + X ) Tr ^ (9JiiJi 2 ■ ■ ■ JJ ■ (24) 

n=0 j 

Generators of the Lie algebra form a basis of, say, left invariant vector fields, so that the 
above formula is really a way of computing Lie derivatives with matrices (along the vector 
field X\ g = L g *X). 

Using the above formulae, we can give a meaning to the derivatives of the propagator. 
To keep notations simple at this stage (the details are given in the next sections), we write 
its expansion in a very symbolic form: 

00 

C({gi, 9 2,g3h ih,h e x \ h e X3 }) ^J2^ X T C {n) ({91,92,93}; , (25) 

n=0 



We will systematically use the sum over repeated indices for contractions on the Lie algebra. 
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where (X) n stands for all products of order n of the form: fc , ^_ k y X % 2 ■ ■ -X^X^ 1 ■ ■ ■ X% n h , 

and the corresponding derivatives evaluated at X2 = X 3 = 0. 

Notice that the apparent symmetry breaking between the three strands initially intro- 
duced appears here as there are no derivatives of the propagator with respect to variables of 
the strand 1. 

The form of the expansion - Finally, the localization of the amplitude ( 117]) leads to: 



/oo 
d 3 X 2 d*X 3 [£ ~ {XT C^({g s }; {g s }) 



x 2 +x 2 +(x 2 -x 3 ) 2 



n=0 

S>3 (X2 ,-^3) 

(i{X 2 ,X 3 )e-~t . (26) 



Usually, the saddle point method is applied to evaluate a quantity as an asymptotic series, 
here in powers of A t . In that case, at each order we would encounter a lot of terms, 
mixing the expansions of the measure /1, of the remainder S> 3 , and that of the propagator. 
However, this is not what we are looking for. Instead, we will organize the series according 
to the order of the derivatives of the propagator. This means that we push the sum over n 
in the above formula out of the integral, together with the derivatives of the propagator: 

Ag;t({9sh{~9s}) * J2 C{n) ^i,92,93h{~9i,h,9s}) pf ] , (27) 

n=0 

with p^ being: 

/ \ C „ „ x 2 +x 2 +(x 9 -x 3 ) 2 s >3 (x 2 ,x 3 ) 

pS n) = A? J d 3 X 2 d s X 3 (X) n e- 3 i »(X 2 , X 3 )e-^~^ . (28) 

It is this coefficient that we want to expand into powers of A t . Actually, our goal is more 
humble than that, since we just want to know whether the leading order of pf^ converges or 
diverges with A t ! 

The components of pf^ satisfy some rotation invariance that we now describe. Let Hi 
denote the Lie algebra su(2) seen as the vector space R 3 equipped with the standard vector 
representation of SU(2). Let us also make explicit the vector indices of p\ n ^: 

I \ ■ f „ „ ■ • X 2 +X 2 + (X 9 -X 3 ) 2 s >3 (x 2 ,x 3 ) 

p (n)n.., n = A? / tf X% #X*X*-.-X* e- j p(X 2 ,X 3 )e-^-^- , (29) 



where each insertion X % is either X\ or X\. It turns out that p^> %1 --- %n j s a invariant tensor 
in "Hf n , that is an intertwiner from "Hf n to C 

Lemma 1. Let g G SU(2) and R{g) denote its matrix in the vector representation. Then, 

n i (^)---^ n (^)A (n)il '-* = P t ( " )il --- i "- (30) 

The proof is straightforward. First, we know that the Haar measure dk 2 dk 3 on SU(2) 2 
is invariant under the conjugation of k 2 and k 3 by g. This implies that the measure 
d 3 X 2 d 3 X 3 p(X 2 , X 3 ) is invariant under the rotation by R(g). Then, we apply the same 
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reasoning with the small t heat kernel approximation. The distances on the group: l^l, \k 3 1 
and jA^A;^ 1 ] are invariant under conjugation by g, and thus their expansion around the saddle 
point (the right hand side of (1201) ) is invariant under the rotation by R(g). Thus, a change 
of variables (X 2 ,X 3 ) H- (R(g)X 2 , R(g)X 3 ) in the above definition of p^ n --- %n l ea ds to the 
desired result. 

The zeroth order: mass renormalization - The leading term to the zeroth order 
n = is simply obtained by setting /ie _5 ^ 3//4t ~ 1, so that: 

p (Q) = A 9- 6 ^(o) + o(A?) (31) 

where the constant is basically the inverse square root of the determinant of the Hessian, 
and is independent of t. Since the first correction to the propagator coming from the graph of 
Figure [3] is a divergent factor A 3 times the bare propagator, it gives the mass renormalization: 

A g;t ({g s };{g s })~A 3 t KM C({g 1 ,g 2 ,g 3 };{g 1 ,g2,93}) • (32) 

Let us go now to the derivatives of the bare propagator ! 

The first order is vanishing - This is a simple corollary of Lemma HJ There is a single 
invariant vector in the zero vector. Hence, 

P? = 0, (33) 

to all orders in A t . 

Even orders and the need for a wave function renormalization - Since the even 
Gaussian moments are non-zero, it is clear that the leading term to p[ 2m ^ corresponds to 
approximating the measure to fi ~ 1 and neglecting the expansion S> 3 . Further, the (2m)- 
th Gaussian moment picks up an extra factor A^ 2m compared with the above zeroth order: 

(am) ^ A 3- 2m J_ f ^X 2 d 3 X 3 (X) 2m e -(x!+x! + (x 2 -x :i f) (34) 

7T m J 

Let us first focus on the scaling properties. It turns out that A^~ 2m goes to zero with t — > 0, 
as soon as m > 2. However, for m — 1, the 2-point function receives a divergent contribution, 
coming with second derivatives of the bare propagator: 

Ag;t({9sY,{9s}) * C({g s };{g s }) (a^°) +o(A f 3 )) +C^({g s }; {~g s }) +o(A t )) . (35) 

This is one of the main result of our analysis, and the details of the terms which are contained 
m are reported to the following section. 

Odd orders are finite - Finally, we can show that there are no divergences in front 
of odd derivatives of the propagator. Since odd Gaussian moments are zero, non-zero con- 
tributions to pp m+3 ' ) involve the expansions of the measure /i and 5*>3. More precisely, we 
have to consider odd terms in their expansion, to create an even term with (X) 2m+3 . A key 
observation is that the expansion of the measure does not contain linear terms in X 2 ,X 3 , 

V(X 2: X 3 ) = 1 + o(X 2 ) + o(X 3 ) . (36) 



10 



Then, the expansion of S> 3 begins like: 

e" £M ^ = l-^ £ S§irxiXlX5 + o{X°). (37) 

s,u,v=2,3 

Thus, the first non-zero terms in pf' comes from the contraction of S^ uv with Gaussian 
moments of order 6. The latter scale like AjT 3 , so that the divergences are fully exhausted 
through (1351) . 

3.3 The second derivatives of the propagator 

Since the first radiative correction exhibits a divergence which cannot be reabsorbed into the 
bare parameters of the action, the corresponding terms will need to be added to the action 
from the beginning (see the discussion in Section [H]). So it is worth giving some details on 
these terms which involve the second derivatives of the propagator. In particular, it is not 
clear from our saddle point analysis why the corrections should be symmetric in the three 
strands, and not only the strands 2, 3. 

(2) 

Symmetries of the Gaussian moments - The coefficient p\ in front of the second 
derivatives of C contains three types of terms: 

p™ = A t i- J d 3 X 2 d 3 X 3 XI XI e -<*? + ^+(* a -*.) a ) , for 8 = 2, 3, (38) 
PiS =At lj ^ ^ X 2 X( e -(*M+(*.-*0 a ) . (39) 

The indices i,j are indices of 3- vectors (thus transforming under the vector representation of 
SU(2)), and are chosen for convenience to be those of an orthonormal basis (the covariance 
can be restored using the metric instead of Kronecker delta in the following). 

Let g G SU(2), and denote R(g) the matrix of the vector representation. Lemma [T] states: 

m\m'i$? =p&s ■ (40) 

Thus, integrating over all rotations, we see that ^ is zero unless i = j. Furthermore, the 
three diagonal terms are equal: 

Pt,rs - ° Q Pt,rs ■ l 4i J 

Using the symmetry under the exchange of X 2 ,X 3 , one also gets pf^2 = pffis- 

The derivatives of the propagator - Now let us come to the consequences of these 
symmetries in the form of the radiative corrections. More precisely, they teach us that 
one should not consider arbitrary combinations of the second derivatives of the propagator. 
Indeed, let us focus on the case r = s = 2 first. Since the Gaussian moments produce a 
Kronecker 8 % i , we consider: 

C% = j dh5( gi h~ gi ) 5(g 3 hg 3 ) ]T(2j2 + 1) Tr J2 (g 2 h J, J, & ~g 2 ) 

= J dh8(g 1 hg 1 ) 5(g 3 hg 3 ) ^(2j 2 + 1) [-j 2 (j 2 + 1)] Tr i2 (g 2 hg 2 ) (42) 



.12 
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in which we have recognized the Casimir J 2 = —j(j + 1) in the representation of spin j. A 
similar formula holds for the strand 3. 

We also have to consider crossed terms in with one derivative on the strand 2 and 
one on 3, 

Cl%= Yl (2ji + l)(2j 2 + l)(2j 3 + l) J dh Tr^hh) [s ij Tr i2 (g 2 h J { g 2 ) Tr h {g 3 h J 3 g 3 ) 

31, 02,33 

(43) 

It turns out that the integral over h here is crucial. Indeed it projects onto the sector which 
is invariant under the group action, so that the standard recoupling theory of SU(2) can 
be applied The insertion of <5 iJ Jj <S> Jj between the strands 2 and 3 can be seen as a 
grasping in the vector representation (which is obviously the spin 1). This grasping actually 
has a diagonal action (when the integral over h is really taken into consideration), 



I 32 


J2 


1} 


f dh 


[33 


33 







= E l[Mh + ^+ 3s(j3 + l) -3i(3i + I)] fdh J] (2j s + l)Tr js (g s hg s ).(U) 

31,32,33 s=l,2,3 

In the first line, we have introduced some coefficients Nj = \J (2j + + 1) and a Wigner 
6j-symbol with a spin being 1. The latter is then explicitly evaluated to arrive at the second 
line. 

The final evaluation - So the result of this operator is actually a linear combination 
of Casimirs, which could just be reabsorbed into Laplace operators. But we want to see 
what combination of Laplace operators comes out, and how the symmetry between the three 
strands is restored. 

So we cannot really go further just using the symmetries, and we will make use of the 
explicit Hessian matrix instead. The integral (I38p is indeed really standard and produces 
Gaussian moments given by the matrix elements of the inverse of Hessian (|22|) . So we 
contract the derivatives with the Gaussian moments to get: 

ifHess- 1 )^ + ^(Hess- 1 )^ + (Hess" 1 )^^ (45) 

= {-^(32 + 1) - 33(33 + 1) + 2 [32(32 + 1) + 33(33 + 1) - ji(h + 1)] } 

(46) 

■ dh H (2j s + l)Tr js (g s hg s ) 

3 s=l,2,3 



31,32,33 



= ~ E [3i(3i + l) +32(32 + 1) +33(33 + 1)] I [dh H (2j s + l)Tr js (g 8 hg s ) . (47) 

31,32,33 s=l,2,3 

Thus, the explicit form of the Hessian has restored the symmetry between the three strands. 

Moreover, the above result can be recast in terms of a differential (Laplace) operator act- 
ing on the bare propagator. Indeed, it is known that the characters are eigenfunctions of the 
Laplace operator A on the 3-sphere, with eigenvalues being the Casimir of the representation, 

A ^(g) = -j(j + 1) Tr^g) . (48) 
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Thus, we have our final formula for the expansion of the 2-point graph: 



AgMs.Y, {«.})= fA? + o(A 3 t )) C({g,}; {§,}) + A, AT< 2 > 



E * 



w c({&}; 

s=l,2,3 

+ convergent terms . (49) 



This study of the scaling behavior of the 2-point function is clearly suggestive about 
which kind of terms have to be added in the Lagrangian at the very beginning, before re- 
discussing the renormalizability of the 3d model. A Laplace operator on the group needs be 
considered in the dynamics in order to achieve a consistent renormalizability of the Boulatov 
model. Further comments on such a theory are reported in Section [61 

We will find quite generally (in any dimensions greater than 3 see the next section) in 
Section [3] that there are divergences with second derivatives of the propagator, which cannot 
be reabsorbed with a counter-term in the initial action. But we have not worked out their 
precise form like the 3d Boulatov case. This is work in progress. 



3.4 /^-dimensional extension 

The similar graph Q (see Fig. [TJ in dimension D has the following amplitude: 



M{9sY,{9s}) 



I) 



\dh e 



D 



\s (gMgs)- 1 ) 



s=l 



n 6 (k v 1 ) 



(50) 



l<a<b<D 



where {g s }f=i and {g s }f=i are the external group elements. 

The computation follows exactly the steps of the previous case. One can reabsorb hi into 
each h s , s — 2, . . . , D, so that it disappears from all internal faces, but appears instead in all 
external faces. The corresponding integral is then seen as the group averaging defining the 
bare propagator. The regularization proceeds by using the small time behaviour of the heat 
kernel on all internal faces. Then, one perform a saddle point approximation around h s — 1, 
for s — 2, . . . , D. One expands h s = e x " into powers of X s e su(2), and the propagator is 
correspondingly expanded. This yields: 



A S A{9sh{~9s}) = J2 C{n) ({9sh{~9s}) 



Pi 



(n) 



(51) 



n>0 



with 



Pi 



(n) 



A 



\D(D-1) 



D 



\d 3 X s 



{X) n exp 



, s>2 



E 

2<r<s<D 



(X r — X s ) 2 R{X S 



(52) 

Here R(X S ) is everything which comes from higher order expansions of the Riemannian 
distances on the manifold and of the Haar measure. Obviously, the leading order of each 
pf^ is evaluated with R(X S ) ~ 1. 

(n) 

Again, the behaviour of p\ in terms of its vector indices is easy to obtain using the 
rotation symmetry of the integrand. The result is that it is an intertwiner, from Hf 11 to C 
The expansion then produces: 
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At the zeroth order 



AgA{9s}]{g s })=A 



§(D-l)(D-2) 



C({g s };{g s }) ff<°> + o(A| 



|(D-l)(D-2) 



) 



(53) 



Not so surprisingly, this is finite for D = 1 or D = 2. 

• The only rotation invariant vector in Hi is zero. Hence, pf^ = 0, to all orders in A t . 

• For the second derivatives of the propagator, it is found again that as a tensor in the 
vector indices, is proportional to the metric, say 5 % i in an orthonormal basis of 
Hi. The Gaussian moment brings a factor A;" 2 , so that: 



The Laplace operators come from the diagonal insertions, X\X 3 S in p\ , and the crossed 
derivatives from X % r X{ with r ^ s. 

Beyond this second order, one can unravel subleading divergences, in front of derivatives of 
the bare propagator to even and odd orders, as long as the involved Gaussian moments are 
of order 2k, with: §(£> - 1)(D - 2) - 2k > 0. 

Exactly like in the 3d Boulatov model, a divergence with second order derivatives is found, 
which cannot be balanced with a counter-term of the bare action. Thus, the latter has to 
be amended to include the necessary quadratic operators from the beginning. However, 
unlike in the 3d case, we only know the expression (1541) for these operators, which is not 
explicitly symmetric in the exchange of the D strands. We nevertheless know that it is 
symmetric, though not in an obvious way. The symmetric expression in 3d was obtained by 
explicitly computing the action of the graspings <9( r ) ® d< s \, using the rotation invariance of 
the propagator (see equation ( l4*4"j) ). Preliminary computations in the four-dimensional case 
exhibits a sum of the D Laplace operators. But it is also possible that in D > 4, there is 
some room for other rotation invariant operators. 

4 Arbitrary 2-point graphs 

We would like to generalize the previous method to any graph of the 2-point function, that 
is with an arbitrary number of internal vertices, and to arbitrary dimensions. The recipe we 
used can be summarized as follows: 

• Use a gauge fixing procedure to reduce the number of integrals. 

• Localize the integral thanks to the (regularized) delta functions. 

• Evaluate the external propagators and their derivatives on the saddle point. 



AgA{g s h{g s }) = A 




+ A ^ (2) (^E A ^+ E dlr)®d(s)*)c({g s };{~g s }) + o(At 2 )\ . (54) 



s=2 2<r<s<D 
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Since the amplitude is again of the form of (jTOj) . these three steps can actually be carried 
out for any graphs of the theory. 

However, it will generically be impossible to extract the bare propagator and its deriva- 
tives like we did before, for the following reason. The localization, on the second step, takes 
place on the set J 7 of flat connections on the 2-complex. This set (see below) is determined by 
the fundamental group tti{G) and is generically an algebraic variety (when G is an algebraic 
group). Exactly like in the weak coupling limit of 2d Yang-Mills, [25], one can recast the 
amplitude as an integral over the normal bundle to J 7 , then perform the integrals over the 
normal fibers (the transverse directions to J 7 ) thanks to the Gaussian behaviour of the heat 
kernels in these directions, [T6] . 

The net result is that when applying the third step in the above recipe, the evaluation 
of the external propagators does not produce the bare propagator of the model, but instead 
one faces some averaging of it over the set J 7 of flat connections on T. This set is potentially 
non-trivial, which means that the group elements on the edges of T to be integrated cannot 
be localized on the identity. This is obviously in contrast to what happens in the previous 
section with the presence of a single saddle point, k 2 = = t. 

Thus, to get the same kind of expansion as that of the previous section, we will later 
restrict attention to graphs with a trivial internal fundamental group, that we call simply 
connected 2-point graphs. 



4.1 The gauge- fixing procedure 

To understand the gauge fixing procedure, it is useful to remember how the amplitude ({TO]) 
is built as the partition function (with external data) for a system of connections on a 2- 
complex. Let L(Q) be the set of internal lines (edges). A discrete connection A is a map 
from L(Q) to G, or equivalently, a collection of group elements assigned to the lines, 

A = (hi)i eL (g) ■ (55) 

These group elements can be thought of as parallel transport operators, or holonomies, along 
each line, between their source and target vertices. The curvature is expected to describe 
the effect of parallel transport along a closed path in Q. So if F mt (Q) is the set of internal 
faces, we consider the oriented product of the holonomies along the boundary of each face. 
This defines a map: 

H : A^(H f =Y[h; if ) feFint{gy (56) 
ledf 

where = ±1 depends on the relative orientation of the line / and the face /. 

We assume that the three open, external faces are each identified by a strand s = 1, 2, 3, 
which goes from one external leg to the otheJ§. The parallel transport operators along them 
are inserted between the boundary variables, say g s on the left and g s on the right, 

H 8 {g„g„A)=g a [l[ hf f ] . (57) 

iedf s 



6 This is the hypothesis of our paper which is not generic in ordinary group field theory. Still, it docs hold 
without restriction in the colored model [2T1. 
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The amplitude then reads: 



Ag({gs};{9s}) = 




Like in the previous section, we regularize the internal delta functions with heat kernels that 
are approximated by their small time behaviour: 

Ag;t({9sY, {~9s}) = Af itnG)|i?iat| I X\dh J] e JJh ^ JJ S^H s (g s ,g s ,A)^ , (59) 

1&L(G) /GF int (0) s=l,2,3 

where F mt is the number of internal faces. 

We now come to the gauge fixing. It seems at first that the standard procedure is not 
available due to the external legs, since there are boundary variables {g s ,g s ) which are not 
integrated. However the integrals over the group elements hi on the links adjacent to the 
external vertices are sufficient to get the standard gauge invariance. Let us say: 

Lemma 2. Inserting and integrating over some parallel transport group elements h, h on the 
external legs of any 2-point graph leaves the amplitude invariant: 

[ dhdhA g . t {{g s h};{g s h})=A g . t {{g s y,{g s }). (60) 

(The result actually holds even without integrating). 

The proof is straightforward. Consider one of the external vertices, say that on the 
left u cxt . Assume without loss of generality that the three internal adjacent links to t> cxt 
are oriented outwards. Let h G G. From a connection A, define a new connection 
by changing the group elements on the three adjacent links to w ext to ki = h[h~ l . Let 
us investigate how it modifies flB^l) . This change of variables leaves the measure invariant: 
[Jj dhi = Ylidki. Quite clearly also, the curvature around the internal faces is unchanged: 
Hf(A) = Hf(^A). However, it does change the holonomies along the external faces to: 

H a (jg 8 ,g 8 WA) = H a (jg 8 h,g 8 ,A). (61) 

Since the amplitude is all in all independent of h, we can integrate it as well, with the 
normalized Haar measure (G is compact). Then, repeat the process on the second external 
vertex {; cxt . 

What we have gained is that the integrand is now explicitly invariant under a natural 
group action acting on all the vertices (and not only internal). This is an action of G^ v ^ 
parametrized by a collection of group elements attached to the vertices, g = (g v ), which 
changes a connection A as well as h, h to: 

{9) A = {g m h ig ^) leL(gy M h = hg -l, ( 9 )~ h = ~ hg -l . (62) 

This symmetry can be used to gauge- fixed the integrand through a standard procedure, |27j . 
which amounts to setting hi — 1 on every line of a maximal tree T on Q. The tree touches 
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all vertices of Q, without forming a loop, and does not contain the external lines. Notice 
that the introduction of the variables h, h in Q6Up is a matter of convenience. That makes it 
possible to apply the result of [27] as it is. Without introducing the additional variables h, h, 
the same result could have been achieved, but with a less straightforward proof, in which 
what happens on the external vertices v ext , v cxt should have been carefully studied. 

Let us write the deformation retract of Q along the tree T . Thus is a 2-complex 
with a single vertex, with (|L(£/)| — | V(C?)| + 1) lines, while the number of faces is unchanged. 
A connection A is redefined to be collection of group elements on Q^> . The result of the 
gauge-fixing process leads to: 



3 

dhdh\{8[g s h[\{ ti} f \{$.h)- 1 



• (63) 



It is convenient to get rid of one of the two integrals, say over h in this formula. This is easily 
achieved by changing the variables according to: hi t— > h^hih, and h i— > hh. The remaining 
integral over h is very useful: it can be viewed as the group averaging which defines the bare 
propagator in (J7J ! Thus, we reach the following 

Lemma 3. Gauge-fixing on a tree of a 2-point graph - The amplitude of a 2-point graph Q 
can be reduced after retraction of a maximal tree to Q^ T \ to: 

A g ,(i9s}; {~g s }) = Af mG)|Fint| / II dh > II ^^C^II^Ka})- (64) 

ieL(g(T)) feF^(g) iedf s 

The map H which sends the elements hi attached to the remaining edges to group elements 
on faces gives a presentation of the fundamental group iri(G) (see below). 



4.2 The generic saddle point analysis 

The equation (164)) gives us the opportunity to get the leading order of the amplitude for an 
arbitrary 2-point graph, using the results of [16]. Indeed, the regularized delta, here with a 
Gaussian behaviour, clearly enforces a localization on the set of flat connections T defined 
as: 

F = H"\t). (65) 

This set admits a nice (and well-known) geometric description in terms of the 2-complex 
Q^ T \ Indeed, one can read a presentation of the fundamental group tti(Q) from it: there 
is one generator per edge of Q^ T \ and one relation per face. The generators are in one- 
to-one correspondence with the group elements to be integrated in (]64p . and the relations 
among the generators are exactly the conditions enforced by the (Gaussian regularized) Dirac 
delta in the amplitude (see details in [16J). So as expected (that was the initial motivation 
of Boulatov [T]), the amplitude is localized around the set J 7 which is determined by the 
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fundamental group. This is the set of homomorphims of the fundamental group into G, also 
known in the mathematical literature as the representation variety of 7Ti(r) into G, 

T = Hom(7n(r),SU(2)) . (66) 

Exactly like in the weak coupling limit of 2d Yang-Mills, [25], one can then recast the 
amplitude as an integral over the normal bundle to J 7 , then perform the integrals over the 
normal fibers (the transverse directions to J 7 ) thanks to the Gaussian behaviour of the heat 
kernels in these directions, [T6] . 

Thus, | Hf (A) | 2 can be linearized in the directions which are transverse to J 7 , around each 
flat connection $ = (0z); e L(g( r ))- Let us denote X$ the variations along these directions. 
Then, the Hessian is given by the quadratic terms in 

\H f (A)\ 2 = \\dHz(Xz)\\ 2 + S> 3; *(X $ ) . (67) 

/ 

Here || • || is the norm coming from the Killing form on g pint . 

To extract the leading order, it is sufficient to evaluate the inserted propagator on J 7 . So 
the integrals over the normal fibers over each flat connection $ can be carried out exactly 
like in [16]. This gives: 

A Q; t({9sh {9s}) * A? {G) I voV($) C({g s J] </>?'}, {~g s }) , (68) 

Jt iedf s 

where the divergence degree^ VL{Q) is: 

fl(g) = (dimG)|F int | - dim(ker dH*) 1 - . (69) 

The volume form vo!jf($) can be expressed using the Reidemeister torsion of the 2-complex 
Q . Generically, while this formula indeed factorizes the leading divergence, it also shows that 
the bare propagator cannot be extracted. The latter is instead somehow averaged over F. 

For instance, if the fundamental group of Y is that of a lens space, vri(r) = Z p , for p e N, 
then there are non-trivial (reducible but not central) flat connections, and the moduli space, 
J-/G, has a finite number of points. The volume of the orbits of the group action can be 
factored out, so that one has: 

■\;A{!,,}: {~9s}) * ^ iS) E "([*]) C({ 9s II Mi'h {~9s}) ■ (70) 

[$}€T/G l&dfs 

While these graph amplitudes do not admit a direct expansion of the kind ([3]), it may be 
still possible to extract some mass and wave-function renormalization using a prescription 
and a scale around which the expansion is enforced, like one usually does in ordinary scalar 
field theorjl^. This certainly deserves to be further investigated. 

7 This formula is an integral over the non-singular subset of J 7 , on which dim(ker dH®) 1 ^ is constant. 

8 In ordinary field theory, a 2-point graph does not take the simple form ([3|). Nevertheless, the mass 
renormalization 5 m is typically extracted by evaluating the amplitude at a given external scale fi, and the 
wave-function counter-term 8z by evaluating the derivatives of the amplitude with respect to the external 
momenta, at the same scale \x. 
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5 Simply connected 2-point graphs 

We will not go further into the study of generic 2-point graphs, and instead focus on a nice 
class of graphs for which it is possible to extract the bare propagator at the leading order, 
and its derivatives at the sub-leading orders: the simply connected 2-point graphs. So from 
now on, we assume that Q^> gives a presentation of the trivial group, 

*i{Q {T) ) = ((a^gmy, (r, = J] «?0 = W ■ ( 71 ) 

ledf 

Thus, there is a single flat connection, namely the trivial one, 

V/ E L(£ (r) ) hi = l. (72) 

Combining this fact with the expression (I64p of the amplitude allows to apply the specific 
result of [IT] concerning simply connected graphs. We briefly outline the reasoning. We 
will perform a saddle point approximation of (151)) . The Hessian is generically given in ( 167)) . 
But on the trivial connection, the differential of the curvature map H simplifies, and is 
actually the first coboundary operator of the 2-complex (with the external edges and 
faces removed): 

dHx = 5 1 (g {T \R) . (73) 
The non- degeneracy of the Hessian is obtained through: 

ker dH x = H 1 (£ (r) , R) = {0} , (74) 

where the first equality follows from the fact that has a single vertex, and the second 
from the Hurewicz theorem. 

Then, the following steps are exactly those we performed on the simple graph in Section 
13.21 The group elements are expanded into powers of Lie algebra elements, hi = e Xl , and 
the measure becomes: 

JJ dh t = J] d 3 X L /i({X,}) . (75) 

ieL(g(T)) ieL(g<T)) 

We again expand symbolically the propagator: 

OO 

C({g s J] eeifXl hi9s}) = W n C (n) ({91,92, 93}; {guhrgs}) , (76) 

where (X) n stands for all products of the form: YIi 11^=1 Xf 1 , with Y^iPi = n -> the 
corresponding derivatives, evaluated at Xi = 0. 

So the amplitude takes the form of an expansion into derivatives of the bare propagator, 
as we wanted: 

00 _ 

A g ., t ({9s};{g s }) * J2- Cin) (i9i, 92,93}; {91,-92,-93}) P t ] , (77) 

n=0 

with pf^ 1 being: 

ftW = Aro)l F«, j^ Xi fAT e _M^ M{ lfi})e _i4™. (7g) 

J 1 
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• Using the same method as that of Lemma[T] (taking advantage of the rotation symmetry 
of the integrand), it is found that pi seen as a tensor in "Hf n (in its vector indices) is 
actually an invariant tensor. 

• The divergence degree at the zeroth order is: 

n(0) = (dimC) (\F^{G)\ - \L(G)\ + \V(G)\ - l) (79) 
= (dimG)( X (r)-l) . (80) 

In particular, if T is the 2-skeleton of a cell decomposition of a closed orientable 3- 
manifold, then its Euler characteristic is the number of 3-cells of the cell decomposition, 
and we recover some standard result of the field. The amplitude is thus: 

Ag.,t({g s }; {~9s}) = A? iG) K^ C ({g s }, {g s }) + o(A? {G) ) , (81) 

where the constant comes from the Gaussian integral. 

• There is no term with first derivatives of the propagator, since pf^ = to all orders 
(because the only invariant vector in TLi is the zero vector). 

• Next we look at the leading order of pi . It is obtained by evaluating pe~ s ^ : ^ 4t ~ 1 
on the saddle point. The Gaussian moment of order 2m brings a factor A^~ 2m , so that 
pl™^ behaves like: 

Since, the divergence degree Q(Q) can be arbitrarily high, it also means we can have 
an arbitrarily high number of derivatives of the propagator with divergences. More 
precisely, this happens for all m such that: 

2m < fi(£?) . (83) 



6 Towards a new scale analysis of the model 

We have studied graphs of the 2-point function in the D-dimensional extension of the Boula- 
tov model. In the case of a trivial fundamental group (that we call simply connected graphs), 
we have been able to recast the amplitude as an expansion whose zeroth order term is the 
bare propagator, and higher orders are given by its derivatives. We have been mainly in- 
terested in finding which derivatives come with a divergence (a positive exponent of the 
ultraviolet cut-off A). Like in identically distributed matrix models, we found that arbitrar- 
ily high orders of derivatives of the propagator receive divergent factors, as the considered 
graphs become larger. 

For the first radiative correction to the propagator in D = 3, at the second order in the 
coupling constant, we have found two divergences: one renormalizing the mass, and one in 
factor of second derivatives of the propagator (and derivatives of higher orders in D > 4). 
In ordinary field theory (say the latter renormalizes the operator <frd 2 <f), and is usually 
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reabsorbed by a renormalization of the wave-function. But since the initial action (and hence 
the propagator) does not contain derivatives of the field, the divergence we observe cannot 
be reabsorbed as a counter-term from the action. This shows that a renormalizable Boulatov 
model requires to include such terms in the bare action. 

So to make sense of the model, one could either consider the large N limit proposed in 
[20] . or consider a new action in which the quadratic part in the field is of the form: 

S k in[4>] = / [JI d 9i] $(91,92,93) I ^2 A (s) + m 2 J 0(01,02,03) , (84) 

•' i=l \s=l J 

As shown from the computation 044p . terms which contract a derivative on g\ with another 
on 02 can be reabsorbed into the Laplace operators. Further, it seems natural to choose 
the same coupling for the three strands of the propagator, as it actually happens in the 
first radiative correction, (j4"9j) . Such a theory has already been considered in the literature, 
by Di Mare and Oriti [7] (and more references therein) for instance. Classical solutions, 
with interesting perturbation properties and some effective dynamics have been investigated 
there. This leads to a scenario consistent with the emergence of matter fields as a phase in 
the context. Thus, besides the fact that this is one of the most natural action, and other 
interesting features as we have just mentioned, the term we propose to add appears in the 
present study as an important quantity concerning the renormalization. 

Still, a detailed treatment as a quantum field theory, starting with power-counting results 
like those obtained in the (standard) Boulatov model, are missing. Let us sketch very 
basically the analysis of the propagator. After Fourier transform on SU(2) 3 , the propagator 
is given by: 

C(ji,h,h) = ^ 3 ~, , 2 • ( 85 ) 



J2s=ljs(js + l 



m 



Back in direct space, and after composition with the projector on the sector invariant under 
the diagonal group action, the propagator reads: 



r 3 

C({0 S }, {&}) = dh \{Ph + l)x j A9sh~g: l )]CUi,32, 



Js) ■ (86) 



To prepare a scale analysis, we introduce the Schwinger parametric representation of the 
propagator 



3 



00 

2^ 



C({g s },{~9s})= dhU^ + ^x^hg: 1 )] / df e -*(£.i.W.+D-H» a ) (87) 

J1,32J3 *=1 ° 

POO f 

= / dte~ tm2 / dh T\K t {g a hg- x ) . (88) 
Jo J s=1 



In the second line, we have noticed that the sum over the spin of each strand produces exactly 
the heat kernel on SU(2). Note also that the latter does not appear as a regularization 
anymore. 
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Taking M > 1, we can see that the form of the propagator makes the introduction of 
scales very natural, [2B]. We write C = J2 { Ci, where Cj is the sliced propagator (for the 
massless situation), 



From what we know about the behaviour of the heat kernel, it is clear that at high scales, 
% 3> 1, the time t is small, and hence the propagator non-zero only on a small neighborhood 
of the set (gshg' 1 ) = 1. Taking inspiration in the Euclidean case, we may conjecture the 
following bound for the slice v. 



where A, 5 are some constants. The factor j\/p( 3dimG - 2 ) j s the ultraviolet cut-off, while the 
factor M 2t in the exponential rather plays the role of an infrared regulator. Again, at large 
scales i, the propagator will grow but only if (gghg' 1 ) comes closer to the identity. 

We hope the rather sketchy picture we have just drawn can be further extended to the 
amplitude of the graphs themselves, and that a new relevant "locality principle" can be 
reached which would enable renormalization. 

Finally we note that it is not clear but certainly very interesting to investigate the non 
simply connected graphs. Though they do not directly admit a simple expansion like 03), it 
may be possible to extract some mass and wave-function renormalization using a prescription 
and a external scale around which the expansion can be settled, like in ordinary scalar field 
theory. Mathematically, we know that the insertions of (closed) knots and links in this theory 
produces knot polynomials, so it may also be fruitful to consider insertions of open strands 
like those of the 2-point function. 
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